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Abstract. Global perturbative QCD analyses, based on large data sets from electron-proton 
and hadron collider experiments, provide tight constraints on the parton distribution function 
(PDF) in the proton. The extension of these analyses to nuclear parton distributions (nPDF) 
has attracted much interest in recent years. nPDFs are needed as benchmarks for the 
characterization of hot QCD matter in nucleus- nucleus collisions, and attract further interest 
since they may show novel signatures of non- linear density-dependent QCD evolution. However, 
it is not known from first principles whether the factorization of long-range phenomena into 
process-independent parton distribution, which underlies global PDF extractions for the proton, 
extends to nuclear effects. As a consequence, assessing the reliability of nPDFs for benchmark 
calculations goes beyond testing the numerical accuracy of their extraction and requires 
phenomenological tests of the factorization assumption. Here we argue that a proton-nucleus 
collision program at the LHC would provide a set of measurements allowing for unprecedented 
tests of the factorization assumption underlying global nPDF fits. 



1. Introduction 

Parton distribution functions (PDFs), which define the flux of quarks and gluons (i = q,g) in 
hadrons fi/h{ x i Q 2 )-, P ia Y a central role in the study of high energy collisions involving hadronic 
projectiles. For protons, sets of collinearly factorized universal PDFs are obtained in global 
pQCD analyses based on data from DIS and DY production, as well as W/Z and jet production 
at hadron colliders. These data provide tight constraints on PDFs over logarithmically wide 
ranges in Q 2 and x, and have allowed precision testing of linear pQCD evolution. However, our 
understanding of parton distribution functions in nuclei of nucleon number A, fi/A( x >Q 2 )i is 
much less mature. Knowledge of nuclear parton distribution functions (nPDFs) is important in 
heavy ion collisions at RHIC and at the LHC for a quantitative control of hard processes, which 
are employed as probes of dense QCD matter. Characterizing nPDFs is also of great interest in 
its own right, since the nuclear environment is expected to enhance parton density-dependent 
effects, which can reveal qualitatively novel, non-linear features of QCD evolution. 

Paralleling the determination of proton PDFs, several global QCD analyses of nPDFs have 
been made within the last decade [1-5] based, up until recently, solely on fixed-target nuclear 



DIS and DY data which, compared to the data constraining proton PDFs, are of lower precision 
and cover a much more limited range of Q 2 and x. Such limitation manifests in a particularly 
poor gluon distribution function, since it can only be obtained from logarithmic (5 2 -evolution, 
for which a wide Q 2 -range is mandatory (see [6] for future prospects in DIS on nuclei). 

Moreover, in contrast to the theoretical basis for global analyses of proton PDFs, the 
separability of nuclear effects into process-independent nPDFs and process-dependent but A- 
independent hard processes is not established within the framework of collinear factorized QCD. 

In view of the importance of nPDFs for characterizing benchmark processes in heavy ion 
collisions, it is thus desirable to look for both stringent phenomenological tests of the working 
assumption of global nPDF fits and experimental observables that can be used to constrain 
nPDFs. To this end, we argue that a program of hadron-nucleus collisions at the LHC would 
provide for such tests with unprecedented quality [7]. 



2. Collinearly factorized single inclusive approach to particle production 

To improve on the insufficiency of nuclear DIS and DY data, recent global nPDF analysis [1, 2] 
have included for the first time data from inclusive high-pr hadron production in hadron-nucleus 
scattering measured at RHIC [8-10]. This is the observable we focus on in order to test collinear 
factorization and to determine nPDFs. 

In the factorized QCD ansatz to hadron-nucleus collisions by construction the entire nuclear 
dependence of the cross section for hadro-production resides in the nPDF (see eq.(l) in [7]). It 
is customary to characterize nuclear effects in the parton distribution functions by the ratios 

Rf(x, Q 2 ) = f i/A (x, Q 2 )/f t/p (x, Q 2 ) , (1) 

which show characteristic deviations from unity in global nPDF analysis for all scales of Q 2 tested 
so far. These effects are typically referred to as nuclear shadowing (x < 0.01), anti-shadowing 
(0.01 < x < 0.2), EMC effect (0.2 < x < 0.7) and Fermi motion (x > 0.7). 

Nuclear effects on single inclusive hadron production are typically characterized by the nuclear 
modification factor R^ A 

„h , n daP A ^ h+x L TvA do pp ^ h+x 

<^> = -*^/^^- (2) 

We calculate the single inclusive it production at both RHIC and the LHC within the 
collinear factorization approach. All our calculations use LO PDFs from CTEQ6L [11] with 
nuclear modifications EPS09LO [1] and the KKP fragmentation functions [12]. When the 
collinear factorization ansatz is adopted, the pT-dependence of Rt A traces the x-dependence of 
nPDFs. Although the precise kinematics is complicated by the convolution of the distributions, 
the qualitative dependence of the Bjorken-x of partons inside the nucleus with px of the final 
hadron and rapidity reads x ~ pTe~ y ' / 'y/s. 



2.1. From RHIC to the LHC: testing collinear factorization 

Fig. 1 shows the nuclear modification factor R^ Au for the production of neutral pions in 
- v /snn = 200 GeV deuteron-gold collisions at RHIC, calculated within the factorized ansatz 
at leading order (LO). Results shown in Fig. 1 are also consistent with the NLO-calculation of 
Rj Au in [1]. The RHIC data [8] in Fig. 1 have been used in constraining the nPDF analysis 
EPS09 [1] but they were not employed in a closely related nPDF fit [3], which provides an equally 
satisfactory description of these RHIC data. Therefore, the agreement of data and calculation 
in Fig. 1 is in support of collinear factorization: the enhancement of R^ Au in the region around 
Pt — 4 GeV at mid-rapidity tests momentum fractions in the anti-shadowing region. 



However, qualitatively different explanations of the R^Au measured at RHIC are 
conceivable [13]. While the above calculation accounted for R^AuiV = 0) ^ n terms of a nuclear 
modification of the longitudinal parton momentum distribution only, such nuclear effects can be 
accommodated within a variety of assumptions, namely initial state parton multiple scattering 
(&t broadening) [14] and non-linear low-x QCD evolution [15]. 

Yet, RHIC data being inconclusive about the dynamical explanation of the nuclear 
modification of the spectrum, we argue that repeating the same measurement at the LHC at a 
much higher energy and in the much wider available pt range, will help to disentangle which of 
the theories is correct, or in which kinematical range each one is valid. 





Figure 2. EPS09 LO prediction for the 
nuclear modification factor of neutral pion 
single inclusive production in proton-lead 
collisions at LHC mid-rapidity. All LHC 
rapidities are in the center-of-mass (shift 
to lab Ay ~ 0.47). 
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Figure 1. Nuclear modification factor 
(2) at RHIC mid-rapidity. Data points 
are PHENIX [8]. Red line represents the 
EPS09 LO calculation. The dashed lines 
are the EPS09 uncertainty [1] (black) plus 
data normalization uncertainty (blue). 

As one can see in Fig. 2, the shape of the nuclear modification factor that would be measured 
at LHC mid-rapidity, if the entire nuclear effect in pPb collisions can be factorized into nPDFs, 
is qualitatively different from that observed at RHIC: at more than 40 times higher center of 
mass energy, final state hadrons at the same transverse momentum test O(40) times smaller 
momentum fractions x%. Even tough a shift of the maximum of RpPhi"/] = 0) to such high 
values of pr is a natural consequence in collinear factorization reflecting nuclear modifications 
of PDFs, it cannot be accounted for in terms of A^-broadening, since such models predict a mild 
dependence of the spectrum with the center-of-mass energy [14]. 

When the single-inclusive spectrum is calculated for different values of rapidity the results 
shown in figure Fig. 3 are obtained, to learn that the rapidity dependence of R^ Pb allows one to 
scan the main qualitatively different ranges of standard nPDFs in an unprecedented way. 

The behavior of the nuclear modification factor when saturation is taken into account [15] 
is rather different: a strong suppression factor is found which is incompatible with any of the 
existing nuclear parton distribution sets. Thus, a measurement at the LHC which shows the 
rapidity dependent behavior in Fig. 2 and 3 would be in strong support of collinear factorization 
as opposed to initial state multiple scattering and CGC models. 



2.2. Constraining nPDFs with the LHC 

Despite the perspectives for qualitative tests of collinear factorization, current global analyses of 
nuclear parton distribution functions come with significant uncertainties as we can see in Fig. 4. 
We have compared in Fig. 5 the nuclear modification factor for two nPDF sets, which show 
marked differences: in contrast to EPS09 the gluon distribution of HKN07 [5] does not show an 
anti-shadowing peak but turns for x > 0.2 from suppression to strong enhancement at initial 
scale Q 2 = lGeV 2 . Fig. 5 thus illustrates that within the validity of a collinearly factorized 
approach, LHC data can resolve the qualitative differences between existing nPDF analyses 



10 20 30 40 50 60 70 80 90 

PT (GeV) 



20 30 40 50 60 70 SO 90 

P T (GeV) 



10 20 30 40 50 60 70 80 

Pt (GeV) 



10 20 30 40 50 60 70 80 

Pt (GeV) 



Figure 3. Rapidity dependence of R* Pb (2) for = 8.8 TeV pPb (LHC). The different plots 

scan the dependence from y = —4 (close to Pb projectile rapidity) up to y = 4 (close to proton 
projectile rapidity). Labels indicate whether the nuclear modification originates mainly from 
the shadowing (S), anti-shadowing (AS) or EMC regime. Vertical lines illustrate the rapidity- 
dependent pt range, which can be accessed experimentally with more than iV even t s = 1000 (= 
10) per GeV-bin within one month of running at nominal luminosity. 

and can improve significantly and within a nominally perturbative regime on our knowledge of 
nuclear gluon distribution functions. 
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Figure 4. The nuclear modification factor (1) in 
lead Rf b for three different flavours as given by the 
available nPDF parametrization at LO. 



Figure 5. K£ pb (2) for = 8.8 TeV 

pPb (LHC), from two different sets of 
nPDFs. 
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